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We investigate the eigenvalue spectrum of the staggered Dirac matrix in SU(3) 
gauge theory and in full QCD as well as in quenched U(l) theory on various lattice 
sizes. As a measure of the fluctuation properties of the eigenvalues, we consider 
the nearest-neighbor spacing distribution, P(s). We further study two-color QCD 
at nonzero chemical potential, /i, by constructing the spacing distribution of adja- 
cent eigenvalues in the complex plane. We find that in all regions of their phase 
diagrams, compact lattice gauge theories have bulk spectral correlations given by 
random matrix theory, which is an indication for quantum chaos. In the confine- 
ment phase, the low-lying Dirac spectrum of these quantum field theories is well 
described by random matrix theory, exhibiting universal behavior. 



1 Bulk of the Spectrum 

The properties of the eigenvalues of the Dirac operator are of great interest 
for the universality of certain features of QCD and QED. On the one hand, 
the accumulation of small eigenvalues is, via the Banks-Casher formula, 1 re- 
lated to the spontaneous breaking of chiral symmetry. On the other hand, the 
fluctuation properties of the eigenvalues in the bulk of the spectrum have also 
attracted attention. It was shown in Ref. 2 that on the scale of the mean level 
spacing they are described by random matrix theory (RMT). For example, the 
nearest- neighbor spacing distribution P(s), i.e., the distribution of spacings s 
between adjacent eigenvalues on the unfolded scale, agrees with the Wigner sur- 
mise of RMT. According to the Bohigas-Giannoni-Schmit conjecture, 3 quan- 
tum systems whose classical counterparts are chaotic have a nearest-neighbor 
spacing distribution given by RMT whereas systems whose classical counter- 



parts are intcgrable obey a Poisson distribution, Pp(s) = e~ s . Therefore, the 
specific form of P(s) is often taken as a criterion for the presence or absence 
of "quantum chaos" . 

In RMT, one has to distinguish several universality classes which are de- 
termined by the symmetries of the system. For the case of the QCD Dirac 
operator, this classification was done in Ref. 4. Depending on the number of 
colors and the representation of the quarks, the Dirac operator is described by 
one of the three chiral ensembles of RMT. As far as the fluctuation properties 
in the bulk of the spectrum are concerned, the predictions of the chiral en- 
sembles are identical to those of the ordinary ensembles. 5 In Ref. 2, the Dirac 
matrix was studied for color-SU(2) using both staggered and Wilson fcrmions 
which correspond to the chiral symplectic (chSE) and orthogonal (chOE) en- 
semble, respectively. Here, 6 we study SU(3) with staggered fermions which 
corresponds to the chiral unitary ensemble (chUE). The RMT result for the 
nearest-neighbor spacing distribution can be expressed in terms of so-called 
prolate spheroidal functions, see Ref. 7. A very good approximation to P(s) 
is provided by the Wigner surmise for the unitary ensemble, 

iV(«) = ™s*e-*"* . (1) 

We generated gauge field configurations using the standard Wilson pla- 
quette action for SU(3) with and without dynamical fermions in the Kogut- 
Susskind prescription. We have worked on a 6 x 4 lattice with various values 
of the inverse gauge coupling (3 = 6/g 2 both in the confinement and deconfinc- 
ment phase. We typically produced 10 independent equilibrium configurations 
for each /3. Because of the spectral ergodicity property of RMT one can re- 
place ensemble averages by spectral averages if one is only interested in bulk 
properties. 

The Dirac operator, p = fi + igA, is anti-Hermitian so that the eigenvalues 
X n of ip are real. Because of {^,75} = the nonzero A„ occur in pairs 
of opposite sign. All spectra were checked against the analytical sum rules 
J2 n = and X)a„>o = 3V, where V is the lattice volume. To construct 
the nearest-neighbor spacing distribution from the eigenvalues, one first has to 
"unfold" the spectra. 7 

Figure 1 compares P(s) of full QCD with Nf — 3 flavors and quark mass 
ma = 0.05 to the RMT result. In the confinement as well as in the de- 
confinement phase we observe agreement with RMT up to very high (3 (not 
shown). The observation that P(s) is not influenced by the presence of dy- 
namical quarks is expected from the results of Ref. 5, which apply to the case 
of massless quarks. Our results, and those of Ref. 2, indicate that massive 
dynamical quarks do not affect P(s) either. 



2 



Confinement (3 — 5.2 
ma = 0.05 



Deconfinement (3 = 5.4 
ma = 0.05 




1 2 3 1 2 3 



s s 

Figure 1: Nearest-neighbor spacing distribution P(s) on a 6 3 X 4 lattice in full QCD (his- 
tograms) compared with the random matrix result (solid lines). There are no changes in 
P(s) across the deconfinement phase transition. 

No signs for a transition to Poisson regularity are found. The deconfine- 
ment phase transition does not seem to coincide with a transition in the spacing 
distribution. For very large values of (3 far into the deconfinement region, the 
eigenvalues start to approach the degenerate eigenvalues of the free theory 
given by A 2 = X^=i sia (2Tm li /L ll )/a 2 , where a is the lattice constant, L M is 
the number of lattice sites in the ^(-direction, and = 0, . . . , — 1. In this 
case, the spacing distribution is neither Wigncr nor Poisson. It is possible to 
lift the degeneracies of the free eigenvalues using an asymmetric lattice where 
L x , L y , etc. are relative primes and, for large lattices, the distribution is then 
Poisson, Pp(s) — e~ s , see Fig. 2. 
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Figure 2: Nearest-neighbor spacing distribution P(s) for the free Dirac operator on a 53 X 
47 X 43 X 41 lattice compared with a Poisson distribution, e~ s . 
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Figure 3: Nearest-neighbor spacing distribution P(s) for U(l) gauge theory on an 8 3 X 6 
lattice in the confined phase (left) and in the Coulomb phase (right). The theoretical curves 
are the chUE result, Eq. (1), and the Poisson distribution, Pp(s) = cxp(— s). 



We have also investigated the staggered Dirac spectrum of 4d U(l) gauge 
theory which corresponds to the chUE of RMT but had not been studied before 
in this context. At j3 c « 1.01 U(l) gauge theory undergoes a phase transition 
between a confinement phase with mass gap and monopole excitations for 
[3 < (i c and the Coulomb phase which exhibits a masslcss photon for [3 > 
(3 C . S As for SU(2) and SU(3) gauge groups, we expect the confined phase to 
be described by RMT, whereas free fermions arc known to yield the Poisson 
distribution (see Fig. 2). The question arose whether the Coulomb phase would 
be described by RMT or by the Poisson distribution. 9 The nearest-neighbor 
spacing distributions for an 8 3 x 6 lattice at (3 = 0.9 (confined phase) and at 
(3 = 1.1 (Coulomb phase), averaged over 20 independent configurations, are 
depicted in Fig. 3. Both are well described by the chUE of RMT. 

Physical systems which are described by non-Hermitian operators have 
attracted a lot of attention recently, among others QCD at nonzero chemical 
potential \i. 10 A formulation of the QCD Dirac operator at /i ^ on the lattice 
in the staggered scheme is given by 11 



M Xi y(U,fi) = -^ [Uu{x)f]u{x)5y tX+u - h.c] 

u—x.y.z 

+ Y a [Ui(x) m (x)e^ y ^ - Ujiy^e-^^] (2) 

with the link variables U and the staggered phases r\. For /x > 0, the Dirac 
operator loses its Hermiticity properties so that its eigenvalues become com- 
plex. The aim of the present analysis is to investigate whether non-Hermitian 
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RMT is able to describe the fluctuation properties of the complex eigenvalues 
of the QCD Dirac operator. The eigenvalues are generated on the lattice for 
various values of /i. We apply a two-dimensional unfolding procedure 12 to 
separate the average eigenvalue density from the fluctuations and construct 
the nearest-neighbor spacing distribution, P(s), of adjacent eigenvalues in the 
complex plane. Adjacent eigenvalues are defined to be the pairs for which the 
Euclidean distance in the complex plane is smallest. The data are then com- 
pared to analytical predictions of the Ginibre ensemble 13 of non-Hcrmitian 
RMT, which describes the situation where the real and imaginary parts of the 
strongly correlated eigenvalues have approximately the same average magni- 
tude. In the Ginibre ensemble, the average spectral density is already constant 
inside a circle and zero outside. In this case, unfolding is not necessary, and 
P(s) is given by 14 



Pg(s) = cp(cs) , p(s) = 2s lim 

AT— >oc 



JV-1 



Y[ e„(s 2 ) e 



N-l 

E 



^ n!e„(s 2 ) ' 



(3) 



where e n (x) — X)m=o a;m / m ' anc ^ c = Jo°° ^ s s p( s ) = 1-1429.... For uncorre- 
cted eigenvalues in the complex plane, the Poisson distribution becomes 14 



(4) 



This should not be confused with the Wigner distribution (1). 

We report on simulations done with gauge group SU(2) on a 6 4 lattice 
using (3 — 4/<? 2 = 1.3 in the confinement region for Nf — 2 flavors of staggered 
fcrmions of mass ma = 0.07. For this system the fermion determinant is 
real and lattice simulations become feasible. 15 We sampled 160 independent 
configurations. In the case of color-SU(2), the staggered Dirac operator has 
an extra anti-unitary symmetry 16 and falls in the symmetry class with Dyson 
parameter (3d = 4. 17 However, one can show that the nearest-neighbor spacing 
distribution in the bulk of the spectrum is also given by Eq. (3). 

Our results for P(s) are presented in Fig. 4. As a function of n, we expect 
to find a transition from Wigner to Ginibre behavior in P(s). This was clearly 
seen in color-SU(3) with Nf = 3 flavors and quenched chemical potential, 12 
where differences between both curves are more pronounced. For the symplec- 
tic ensemble of color-SU(2) with staggered fermions, the Wigner and Ginibre 
distributions are very close to each other and thus harder to distinguish. They 
are reproduced by our preliminary data for \x = and \i = 0.4, respectively. 

For [i > 1.0, the lattice results for Pis) deviate substantially from the Gini- 
bre distribution and could be interpreted as Poisson behavior, corresponding 
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Figure 4: Nearest-neighbor spacing distribution in the complex plane for two-color QCD 
with fi in the confinement (left) and dcconfinemcnt (right) phase. The short-dashed curve 
is the Wigncr distribution for the chSE and the solid curve is the Ginibre distribution of 
Eq. (3). 

to uncorrelated eigenvalues. (In the Hcrmitian case at nonzero temperature, 
lattice simulations only show a transition to Poisson behavior for (3 — > oo when 
the physical box size shrinks and the theory becomes free. 6 ) A plausible expla- 
nation of the transition to Poisson behavior is provided by the following two 
(related) observations. First, for large p the terms containing in Eq. (2) 
dominate the Dirac matrix, giving rise to uncorrelated eigenvalues. Second, 
for large p the fermion density on the finite lattice reaches saturation due to 
the limited box size and the Pauli exclusion principle. 



2 Low-lying Spectrum 

We have continued our investigations with a study of the distribution of the 
small eigenvalues in the confined phase. The Banks-Casher formula 1 relates 
the Dirac eigenvalue density p(X) at A = to the chiral condensate, S = 
\{i>i^)\ = lim £ ^o limy^oo np(e)/V. The microscopic spectral density, p s (z) = 
limy^oo p (z/VT,) /V£, should be given by the appropriate result of RMT 18 , 
which also generates the Leutwyler-Smilga sum rules. 19 

To study the smallest eigenvalues, spectral averaging is not possible, and 
one has to produce large numbers of configurations. We present results from 
Ref. 20 for SU(2) theory and the staggered Dirac operator. Both the distribu- 
tion P(A m in) of the smallest eigenvalue and the microscopic spectral density 
p s (z) agree with the RMT predictions of the chSE for topological charge v = 0, 



6 




0.001 0.002 0.003 



0.6 

P s iz) 
0.-1 

0.2 



1416 conf. _ 




13 = 2.0, |y|=10 4 

I i i i i I i i i i I i i 
10 15 20 25 



Figure 5: Distribution P(A m i n ) (left) and microscopic spectral density p s {z) (right) from 
Ref. 20 for SU(2) gauge theory on a 10 4 lattice in comparison with the predictions of the 
chSE of RMT (dashed lines). 



as depicted in Fig. 5. 

Our analog results for U(l) theory are for (3 — 0.9 in the confined phase 
with 10000 configurations on a 6 4 lattice. The left plot in Fig. 6 exhibits 
the distribution P(A m i n ) of the smallest eigenvalue in comparison with the 
prediction of the (quenched) chUE of RMT for topological charge v = 0, 
-P(A m in) = (VE) 2 (A m i n /2) exp(-(V£A m in) 2 /4). For the chiral condensate we 
obtain S ks 0.35 by extrapolating the histogram for p(X) to A = and using 
the Banks-Casher relation. In the right plot of Fig. 6 the same comparison 
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Figure 6: Distribution P(A m i n ) (left) and microscopic spectral density p s (z) (right) from 
our 6 4 lattice data of U(l) gauge theory in comparison with the predictions of the chUE of 
RMT (dashed lines). 
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Figure 7: Density p(A) of small eigenvalues for U(l) on an 8 3 X 6 lattice at (3 = 0.9 (left) 
and at (3 = 1.1 (right). A nonzero chiral condensate is supported in the confinement phase 
of U(l) gauge theory. 

with RMT is done for the microscopic spectral density p s (z) up to z = 10, 
and the agreement is again quite satisfactory. Here, the analytical RMT result 
for the (quenched) chUE and v — is given by 18 p s (z) = z [Jq(z) + Jf(z)]/2, 
where J denotes the Bessel function. 

In U(l) theory, the quasi-zero modes which are responsible for the chiral 
condensate £ w 0.35 build up when we cross from the Coulomb into the 
confined phase. For our 8 3 x 6 lattice, Fig. 7 compares on identical scales 
densities of the small eigenvalues at (3 = 0.9 and at (3 = 1.1, averaged over 20 
configurations. The quasi-zero modes in the left plot are related to the nonzero 
chiral condensate, whereas no such quasi-zero modes are found in the Coulomb 
phase. It would be worthwhile to understand the physical origin of the U(l) 
quasi-zero modes in more detail. For 4d SU(2) and SU(3) gauge theories a 
general interpretation is to link them, and hence the chiral condensate, to the 
existence of instantons. In 4d U(l) gauge theory it has been shown numerically 
that the chiral condensate is mainly due to monopole configurations. 21 

For two-color QCD, numerical simulations of the full theory with chemical 
potential in Ref. 15 have exhibited a chiral phase transition at p c « 0.3 where 
the chiral condensate as the associated order parameter vanishes. In Fig. 8 
we compare the densities of the small eigenvalues at p = and at p = 0.1 to 
0.4 on our 6 4 lattice, averaged over 160 configurations. Since the eigenvalues 
move into the complex plane for p > 0, a band of width e = 0.015 parallel to 
the imaginary axis is considered to construct p(y), i.e. p(y) = dx p(x,y), 
where p(x, y) is the density of the complex eigenvalues x + iy. 

Another interesting question concerns the energy scale E c up to which 
RMT describes the small Dirac eigenvalues in the phase where £ > 0. In 
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Figure 8: Density p(y) of small eigenvalues (see text) for two-color QCD on a 6 4 lattice at 
fi = (left) and at fi = 0.1 to 0.4 (right). The vanishing of the chiral condensate in the 
matter phase, found in Ref. 15 at fi c 0.3, is accompanied by a loss of quasi-zero modes. 

disordered mesoscopic systems, a similar scale is called the Thouless energy. 
The theoretical prediction for QCD is E c ~ f 2 /Y,L 2 (see Ref. 22) with the 
pion decay constant f v , where we have assumed that the spatial extent L s of 
the lattice is not smaller than the temporal extent L t . In units of the mean 
level spacing A = tt/VY, at the origin, this becomes 

u c = f ~ ~flL s L t . (5) 

A convenient quantity from which u c can be extracted is the disconnected 
scalar susceptibility 

x^M = 1 ( E {lXk + m ]( t x l + m) ) A - Jr (g Jx^) a ■ (6) 

The corresponding RMT result for the quenched chSE with v — is given 
by 23 Xrmt = 4u 2 tids S 2 K Q (2su) ^dtI (2stu){ S (l-t 2 )+AK Q (2u)[I Q (2su)+ 
tl (2stu)} - 8stI Q (2stu)K (2su)} - 4u 2 'K 2 {2u)[^ ds I (2su)} 2 , and for the 
quenched chUE with v = we have 24 Xrmt = u 2 [K 2 (u) - K^(u)][I^(u) - 
I 2 (u)]. Here, u = mVY,, and I and K are modified Bessel functions. In Figs. 9 
and 10 we have plotted the ratio 23 

ratio = { X fZ ~ Xr^t) /A (7) 

versus u and u/(L s L t ), for color-SU(2) with staggered Dirac operator at (3 = 
2.0 from Ref. 23 and for the U(l) data computed at /3 = 0.9, respectively. This 
ratio should deviate from zero above the Thouless scale. The expected scaling 
of the Thouless energy with L s Lt is confirmed. 
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Figure 9: The ratio of Eq. (7) for SU(2) gauge theory plotted versus u and u/(L B Lt), 
respectively (in this case, L s = Lt = L). The left plot is for L = 10. In the right plot, the 
data for different L fall on the same curve, confirming the expected scaling of the Thouless 
energy according to Eq. (5). The deviations of the ratio from zero for very small values of u 
are well-understood artifacts of the finite lattice size and of finite statistics. 23 
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Figure 10: Same as Fig. 9 but for U(l) gauge theory. 
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3 Conclusions 



The aim of this paper was to work out two different types of universalities 
inherent in quantum field theories and their interpretation in terms of RMT. 
The first type concerns the bulk of the spectrum of the Dirac operator. The 
nearest- neighbor spacing distribution P(s) agrees with the RMT result in both 
the confinement and the deconfincmcnt phase of pure gauge theory and of full 
QCD, except for extremely large values of (3 where the eigenvalues are known 
analytically. The nearest-neighbor spacing distribution of 4d U(l) quenched 
lattice gauge theory is described by the chUE of RMT in both the confinement 
and the Coulomb phase. Even in the Coulomb as well as in the deconfine- 
ment phase, gauge fields retain a considerable degree of randomness, which 
apparently gives rise to quantum chaos in these theories. 

A general unfolding procedure for the spectra of non-Hcrmitian opera- 
tors was applied to the lattice Dirac operator for two-color QCD at nonzero 
chemical potential. Agreement of the nearest-neighbor spacing distribution in 
the complex plane with predictions of the Ginibre ensemble of non-Hcrmitian 
RMT was found around fi = 0.4. The changes for larger values of fi toward 
a Poisson distribution arc understood formally. The physical interpretation 
requires a better understanding of QCD at nonzero density. 

The second type of universality concerns the low-lying spectra of the Dirac 
operators of both QCD and QED. In all cases considered, one finds that in 
the phase in which chiral symmetry is spontaneously broken the distribution 
-P(A m i n ) and the microscopic spectral density p s {z) are described by RMT. 
The Thouless energy scales with the lattice size as expected. 

In summary, both the bulk and the low-lying Dirac spectrum are governed 
by RMT and the related symmetries. 
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